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Abstract 


In this paper the general problem of a rectangular plate clamped along 
two parallel sides and containing a crack parallel to the clamps is consid- 
ered. The problem is formulated in terms of a system of singular integral 
equations and the asymptotic behavior of the stress state near the corners 
is investigated. Numerical examples are considered for a clamped plate 
without a crack and with a centrally located crack, and the stress intensity 
factors and the stresses along the clamps are calculated. 


1. Introduction 

A long strip containing a central or an edge crack and a rectangular 
block with an edge crack are two of the more widely studied geometries in 
fracture mechanics. The former -generally simulates a single-edge-notched 
or a center-notched specimen and the latter the compact tension specimen. 
These specimens are often used for fatigue crack growth and fracture charac- 
terization of engineering materials. Typical studies using a variety of 
methods such as finite elements, Wiener-Hopf, weight function, Laurent 
series, conformal mapping, boundary collocation, and integral equations may 
be found in references [1]-[13]. In all these and similar studies the 
external loads are usually assumed to be such that the boundary conditions 
may be prescribed in terms of tractions only. However, in some cases the 
loads are applied to the specimen through "grips". In these problems it 
would be more appropriate to prescribe the boundary conditions along the 
grips in terms of displacements rather than tractions. The problem of a 
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rectangular strip without a crack loaded through two rigid grips was con- 
sidered in [14], where it was shown that the contact stresses along the 
grips are singular at the corners. In this paper the problem of a strip 
containing a crack is considered by using a technique which is somewhat 
different than that of [14] and which may be used to solve the clamped 
plate problem for an arbitrarily oriented crack. 


2. Formulation of the Problem 

We will first formulate the problem described in Fig. 1 where the 
clamped strip contains three cracks as shown. The "external load" in the 
problem is the relative clamp displacement u Q . One may express the stress 
and displacement components in the strip as follows: 

b k 

6 r * 

c^-jlx.y) = a sij .(x,y) + z G.j k (x,y,s k ,t)f k (t)dt , 

k- 1 


(i > j— x,y , k-1 ,. . . ,6) , 
b k 

6 f K 

u.j (x,y) = u si (x,y) + z V. k (x,y,s k ,t)f k (t)dt 

k-1 _ 


(i) 


(i-(x,y) , k-1 ... . ,6) 


( 2 ) 


where <t s1 . . and u sl - are associated with an infinite strip without any cracks 
and the functions G.j k and represent the stresses and displacements 
in an infinite plate at the point (x,y) due to dislocations f k distributed 
along the cuts (a k <t<b k ). In terms of the crack surface displacements the 
functions f k are defined as follows (Fig. 1): 

f i(x) = [u y (x,s.+0) - u y (x,s 1 -0)] , (i=l,2) 

f 3^> = W [u x (s 3 +0 ’ y) ' u x (s 3'°’ y);l ’ 
f j^ x ^ = ~ix ^ u x^ x,s j +0 ^ " u x^ x,s j' 0 ^ * (j =4 > 5 ) 

f 6 (y) = Jy [u y (s 6 +0,y) - u y (s 6 -0,y)] , (3a-d) 


- 2 - 



where 


s i+3 _s i * a i+3 =a i * b i+3 =b i’ a i <s i <b i* ^ 1 ’ 2 > 3 ) 


(4) 


and the Green's functions G.^ and V- are given in Appendix A. 

The solution which corresponds to an infinite strip may be expressed 
as follows: 

CO 

u sx (x,y) =5^1 C[A 1 +(f+x)fl 2 ]e'“ x + [-A 3 +(f-x)A 4 ]e“ x 

—CO 

CO 

u sy (x,y) = ^ | [(A 1 +A 2 x)e" ctX + (A 3 +A 4 x)e aX ]e" lay da , (5a, b) 


a sxx (x ’ y)= ^T 


{-[a(A 1 +A 2 x)+^ A 2 ]e" aX - [a(A 3 +A 4 x) A 4 ]e aX }e" laX da ! 


a syy( x, y) = irT A 2 +a ^ A l +A 2 x ^ e A 4 + a(A 3 + A 4 x )]e aX }e laX da, 


a sxy (x ’ y > = * 


{-[a(Ai+A 2 x) +-^ 2 ^- A 2 ]e aA +[a(A 3 +A 4 x)- A 4 ]e aA }e iaA dot. 

o 

(6a-c) 


.-ax. 


K-l 


ax, -lax. 


The unknown functions f-j . ,fg and A^(a),...,A 4 (a) are determined from the 
boundary conditions prescribed on the crack surfaces and the plate boundar- 
ies x=0 and y=H. In the actual problem the crack surfaces are free from 
tractions and along the boundaries of the plate we have (Fig. 1) 

u x (0,y) = 0, u y (0,y) = 0, u x (H,y) = u Q , u y (H,y) = 0, -®<y<» . ( 7 ) 


The superposition technique may be used to solve the problem by first 
obtaining the stress state in a clamped plate without any cracks under the 
boundary conditions (7). The equal and opposite of the stresses found in 
this solution may then be used as crack surface tractions in the perturba- 
tion problem. The ten unknowns f^,...,fg, A^...A 4 would then have to be 
determined from 
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( 8 ) 


u x (0,y)=0, u (0,y)=0, u x (H,y)=0, u y (H,y)=0, -~<y<» , 


o yu. 
3-< o 


a yy (x,y j ) = p j (s) = ■ ^TTT ■ (y j* s J> s ' ,x > a J <s< V 1,2) , 

yU 

CT xx^ X 3 ,y ^ = p 3^ = " 7=T1T ’ ( x 3 =s 3» s= y> a 3 <s<b 3) » 


a X y( x »yj) = Pj(s) = 0, (y j .=s j ._ 3 ,s=x,a j ._ 3 <s<b j ._ 3 ,j=4,5) 


a X y(x 3 ,y} = p 6 (s) = 0 , (s=y,a 3 <s<b 3 ) . 


(9) 


By substituting from equations (2), (5) and (A13)-(A20) (of Appendix 
A) into (8), by inverting the Fourier integrals and by evaluating the rela- 
ted infinite integrals, the functions A-j,...,A^ may be expressed in terms 
of integrals of the following form: 


b o 


A, (a) = E 

1 j=i i 


B..j(a,t)fj(t)dt , (i=l , . . . ,4) , 


(10) 




where the functions • depend also on the constants s-j, s 3 » s 3 , H, and k . 
From (10) and (6) the stress components a ... (i,j=x,y) are then obtained 

5 I J 

as follows: 

b k 

6 

0 sij( x,y ) = E c ij k (x,y,s k ,t)f k (t)dt , (i ,j=x,y) , (11) 

k-1 


where the kernels j k are given in Appendix B. 

The integral equations for the unknown functions f^,...,fg are 
obtained by substituting from (1) and (11) into the crack surface boundary 
conditions (9) and may be expressed as 
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( 12 ) 


6 

I 


k=l 


h nk ( s > t ) f k( t ) dt = P n ( s ) 
a k 


(a n < s <b n J n 1 > • • • >6; a n ~a n _ 3 > 
b n =b n-3 5 n=4 ’ 5 ’ 6) 


where the kernels h nk (n,k=l,...,6) are known in terras of and 
(i J=x,y; k=l,...,6) (see equations (1), (11), (9), Appendix A and Appen- 
dix B). For the crack problem shown in Fig. 1 (i.e., for 0<a k <b k <H, 
k=l,2 and s 2 <a 3 <b 3 <s-j ) from appendices A and B it can be shown that the main 
diagonal elements of the kernels h nk which are contributed by G.. - k are 
Cauchy kernels and consequently h nk can be expressed as follows: 

£ 

h nk^ s,t ) tt(1+k) !>? + k nk^ s,t ^ * ( n » k= l » * • • » 6 ) (13) 

where 6 nk =l for n=k, <$ nk =0 for n^k, and k nk , (n,k=l ,6) are known bounded 
functions in the closed intervals, a k £(t,s)£b k , (k=l,...,6). If one wants 
to solve the problem shown in Fig. 1 for any crack geometry provided 
0 <a k <b k<H, k=l,2 and s 2 <a 3 <b 3 <s-| , one could let 

F k (t) 

f (t ) = * , (k=l,...,6) (14) 

/(t-a k )(b k -t) 

and use a Gaussian integration formula to determine the unknown bounded 
functions F k (t), (k=l ,. . . ,6)[15]. From the definition of the unknown func- 
tions f k (t) as given by (3) and from Fig. 1 it is clear that the integral 
equations (12) must be solved under the following singlevaluedness con- 
ditions: 

. b k 

f k (t)dt = 0 , (k=l,...,6) . (15) 

a k 
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3. The Clamped Plate 

From Fig. 1 it may be seen that the problem would reduce to that of a 
clamped plate if one lets a^=0 and b^=H, (i=l,2). In this case techni- 
cally the integral equations (12) and the singlevaluedness conditions (15) 
are still valid and the input functions P n (s) are still given by (9). 
However, the form of the solution is no longer given by (14). To deter- 
mine the fundamental functions of the system of integral equations (12) 
the singular nature of the kernels k nk (s,t), (n,k=l ,...,6) for a k =0, 
b k =H, (k=l,2,4,5) must be investigated. This may be done through the asymp- 
totic analysis of the infinite integrals giving the appropriate functions 
C ijk in A PP endix B which are related to k nk (see, for example, [14] and 
[15]). First we observe that from (1), (11), (9), (12), (13) and Appen- 
dices A and B, the kernels k nk (s,t) way be expressed as follows: 


k nk (s,t) = k^ k (s,t) + k nk^ s,t ^ ’ ( n >k=l,...,6) > 

co 

kj; k (s,t) = j K nk^ s,t,ct ^ da ’ ( n » k=1 »• • • >6) 


(16) 

07 ) 


where k^ are the nonsingular parts of (see equation (1) and Appendix 
A) and kj^ are the kernels contributed by C.^. If we now denote the asymp- 
totic values of for a-*» by K^ koo (s,t,a) , the singular and Fredholm parts 
of the kernels k^ k may be separated as follows: 


k[j k (s,t) = k nk^ s,t ^ + k nk^ s,t ^ » ( n » k=1 

co 

k nk( s,t ) * f K nk»( s ’ t ’ c ') da ’ ( n * k=1 >-- 


k nk ( s .t) 


CO 

f [Kj k (s,t,a)-^Js,t,o)]do 


nk° 


. . . , 6 ) 

> 6 ) , 

» ( n ^ k - l , . . . ,6 ) 


(18) 

(19) 

( 20 ) 


where the kernels k^ k (s,t) (n,k=l,. . .,6) are bounded in their respective 
closed domains. For S 2 <a 2 <b 2 <s-j , after a relatively straightforward analysis 
it may be shown that 

k nk (s ’ t} = k nn (s ’ t)6 nk ’ < n ’ k=1 » • • • > 6 ) » C21 ) 
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k g g(s,t) = 0 » k 66^ s ’^ ^ * 


( 22 ) 


k nn (Sst) = h C{3_k2)+12s & + 4s2 ^ i+t 

- 27 C(3-k 2 )-12(H-s) £ + 4(H-s) 2 ^] ^ . 

(n=l ,2,4,5) . (23) 

Because of (21) and (22) the dominant part of the integral equations 
(12) for n=3 and n=6 has only Cauchy type singularity and the unknown 
functions f g and f g have the form 

F,(t) 

f,(t) = — j~ , (j-3,6) . (24) 

J At-a 3 J(b 3 -t) 

Expressing now the remaining unknowns as 
F k (t) 

f k (t) = , (k=l ,2,4,5), 0<Re(6,Y)<l , (25) 

k (t-a k ) e (b k -t) Y 

through a function-theoretic analysis it may be shown that (see, for exam- 
ple [15]) 

8 = y, 2 kcostt8-(1+ic 2 )+4(8-1 ) 2 = 0 . (26) 

The characteristic equation (26) is identical to that found in [14] and 
in, for example, [16] for a quarter elastic plane fixed along one of its 
straight boundaries. 

4. The Stress Intensity Factors 

In the crack problem shown in Fig. 1, once the functions f^,..,,f g 
giving the crack surface displacement derivatives or the bounded functions 
F.|,...,F g defined in (14) are determined, the Modes I and II stress 
intensity factors at the crack tips may be obtained from the following 
standard relations: 
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k-j (a i ) 


k 2 (a j> 


= 2u I,*, 


Tim ^(t-a.,.) f.(t) 


. 2y 


J 


✓(b i -a 1 )/ l 2 

W 


i!L 

1+K 


_ 2y 


, (i=l >2,3) 


( a j a j-3’ b j b j-3* j=4 ’ 5,6) • 
k l (b i } = ‘ 1+7 ^(b r t) f,.(t) = - ^ ^— L ~ 


t-+b.- 


k 2 ( b 0 .) - 


yr- lim /2TET-i) f,(t) - - yp- 

1+,; t*b, 3 3 1+,c /(b r a,)/2 

J J J 


Ab r a.)/2 


, (i=l ,2,3) , 


2u ‘ J vu j 


( a j =a j-3* b j =b j-3’ J" 4 ’ 5 ’ 6 ) 


(27a-d) 


These relations are based on the definition of the stress intensity factors 
in terms of normal and shear cleavage stresses at the crack tips of the 
form (Fig. 1) 

k,(a,) = lim /2(a,-x) a (x,s,) , ... 

1 _ * * 


k^a-j) = lim ^(a-j-x) ct^Cx.s^ 

x-va -j 


(28a, b) 


In the "clamped plate" problem it -was shown that at the corners y=s 1 - , 
x=0, x=H, (i=l,2) the displacement derivatives have a singularity of the 
order r~ e , where r is a small distance (in x direction) from the corners 
(see eqs. (25) and (26)). This implies that at the corners the stresses 
may also have a singularity of the order r" B , r now being the radial dis- 
tance from the corner. In practical applications one may be particularly 

interested in the behavior of interface stresses a and a, along the 

clamps. We first observe that the stress state at any point in the medium 
is given in terms of the density functions f-j,...,fg and the kernels G-.^ 
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a nd C.^ (see (1) and (11)). Next, we note that the singular behavior of 
the stresses is determined by the dominant parts of the related kernels 
only and these dominant kernels can be separated through a relatively sim- 
ple asymptotic analysis of and given in the Appendices A and 
B. Also, the singular behaviors of the stresses in the medium at all 
four corners have the same form. Thus, after following a procedure simi- 
lar to that described by the equations (16) through (23), for a small 
distance r from the corner x=0, y=s 1 the stresses along the clamp x=0 
may be expressed as 

H H 

a xx (0 ’ s r r) 5 Wrfcy c | Kii^’^M 1 ^ + K i4(y>t)f 4 (t)dt] , 


o 

H 


o 

H 


where 


CT xy(°’ s l" r ) = w(1+k) ^ K 41 ^’ t > f i + 

o o 

(y=s 1 -r) 

I/S „ *\ _ t(t 2 -r 2 ) ^ 1 r(»c 2 -l)t i 2t(t 2 -r 2 )n 

K n (s r r,t) - + 27 lp+fr~ + > 

„s /. „ _ Kr^t 2 ) 1 r(l+K 2 )r , 4t 2 r 1 

K 14 (s r r,t} - -fiz+flf/- “ 2k L t^+r 2 " (t^) zJ ’ 

„s . r(r 2 -t 2 ) 1 r (l+ic 2 )r 4t 2 r i 

^ 4 r s i"r» t / “ - XF5 ir 2 F “ 7Z 1 F+F~ " (F+r 2 ) 2j » 


K^(y,t)f 4 (t)dt] 


(29a, b) 


,/S _ t(t 2 -r 2 ) j. 1 r (l-»c 2 )t , 2t(t 2 -r 2 )! 

K 4 4(s r r,t) - + 27 t t^F" 


(30a-d) 


The kernels given in (30) are "singular" in the sense that they become 
unbounded as t and r approach the end point (t=0, r=0) simultaneously and 
have the order t~^ or r" 1 . To evaluate the stresses given in (29) we 
first note that by substituting f-j(t) and f^(t) from (25), the leading 
terms of the typical integrals in (29) may be expressed as follows (see, 
for example [15]): 
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H 


tf^t) 2 7T F 1 (0) •] 


dt - 2H B 5in n r B ’ 

0 2 


dr 


tf-j(t) 

t 7 +r 2 ~ 


dt = 


ttBF-j ( 0) 
2H 6 sin f 


^ 4 (t) 

t'^+r 2 - 


dt = 


ttF 4 (0) 1 

2H e cos ^ r 6 ’ 


dr 


H f 4 (t) ^ Tt(8+1)F 4 (0) , 

„ tZ+r " dt " ' 2 H b cos 7 
0 2 


( 31 a-d) 


Thus, it is seen that the leading terms in the expressions of the 
interface stresses a xx and a xy around the corner x=0, y=s-j have a singu- 
larity of the form r"&. Similar to (28), defining now the "stress intensity 
factors" by 

k,(0) = lim fL r e a (0,s,-r) , 
r->0 xx i 

k~(0) = lim SZ r B a Y .,(0,s 1 -r) , (32a, b) 

c r-+0 xy i 

and substituting from (30) and (31) into (29), we obtain 

k-,(0) = - f : v — [(icH-2b)cos ^ F 1 (0)-( > d-3-2g)sin & F 4 (0)] , 

►cH p simrB 

k 2 (0) = - [(<-l+2B)sin ^ F 1 (0 )+(k-3+2b)cos ^ F 4 (0)]. (33a, b) 

icH^sinirB 

The stresses elsewhere along the interfaces x=0 and x=H, of course, may be 
evaluated from eqs. (1) and (11) by using the appropriate kernels from 
Appendices A and B. The stress intensity factors in the other corners may 
be obtained by following a procedure similar to that described in (29)- 
(33). 
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5. The Results 


The first example considered is the problem of a clamped plate with- 
out a crack in order to compare the results with that given in [14] where 
a different technique was used to solve the problem. In [14] the measure 
of the external loads was the resultant force rather than the displace- 
ment u Q used in this paper. Hence, the easiest direct comparison may be 
made for the ratio of the stress intensity factors, k 2 (0)/k-|(0) (see (32) 
for the definitions). For the length to width ratio H/2s-j = 2 and 
<=3-4v this comparison is shown in Table 1. 


Table 1. Comparison of ko (0 )/ki (0 ) obtained in this paper 
(second row) ana in [14] (third row), H/2s]=2. 


V 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

-k 2 /k-j 

0 

0.118610 

0.212352 

0.302276 

0.396751 

0.504213 

(-y^) 

(from [14]) 

0 

0.118614 

0.212355 

0.302280 

0.396751 

0.504209 


For various relative plate dimensions and for <=1.8, the normalized 
stress intensity factors defined by (32) are shown in Table 2. In this 
case (26) gives the power of stress singularity as 0=0. 28882. 


Table 2. Normalized stress intensity factors at the corner x=0, 
y=S] of a rectangular clamped plate for <=1.8, plate 
length=H, plate width=2s-] (see insert in Fig. 3) 

ill! S, ^ llU 


H/s-j 

k (0)/( 0 1 1 

if fnwf 0 1 \ i 

-k 2 (0)/k.,(0) 

1 s 1 

1 

1 

2.030255 

-0.613697 

0.302276 

2 

1.024510 

-0.309685 

0.302276 

4 

0.478944 

-0.144773 

0.302276 

6 

0.318621 

-0.096312 

0.302276 

8 

0.240840 

-0.072800 

0.302276 

10 

0.194190 

-0.058699 

0.302276 
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Figures 2 and 3 show the normalized stresses a (0,y) and cr (0,y) 

^ xx xy 

along the clamp for the case of plane strain and for <=1.8. The cor- 
responding Mode II stress intensity factor is shown in Fig. 4. 

The second example considered is a clamped plate containing a cen- 
tral crack (see the insert in Fig. 5). The normalized Mode I stress 
intensity factor calculated at the crack tips y=+c, x=H/2 is shown in 
Table 3 for <=1.8 and <=2.2. For comparison the results obtained in 
[ 6 ] by using the technique of conformal mapping for <= 2.2 are also 
included in the table. The agreement seems to be fairly good for small 
cracks. However, for larger cracks there is some discrepency. For <=2.2 
and H/2s-j=l the normal stress along the clamp is shown in Fig. 5. 


Table 3. Crack tip stress intensity factor k-j ( 83 ) in clamped 
plate containing a symmetrically located central 
crack, H/2s-j=l, c=(b 3 ~a 3)/2 (Fig. 5). 


2 c 

H 


uu n 

k i( a 3 )/(_ ir 


00 

• 

r— 

II 

<= 2.2 

<=2.2 (Ref. 6 ) 

0.1 

3.15135 

2.60029 

2.60000 

0.2 

3.00550 

2.49602 

2.50000 

0.4 

2.78474 

2.17958 

2.17391 

0.5 

2.34322 

2.00829 

2.01724 

0.667 

1.98095 

1.73227 

1.78663 

0.8 

1.71592 

1.52883 

1.64815 


REFERENCES 

1. W.K. Wilson, in Methods of Analysis and Solutions of Crack Problems, 
ed. G.C. Sih, Noordhoff, Leyden (1973) 484-515, 

2. W.T. Koiter, Journal of Applied Mechanics Transactions 32, ASME 

(1965) 237-242. 

3. H.F. Bueckner, Zeitschrift fur angewandte Mathematik und Mechanik 51 
(1971) 97-109. 

4. M. Isida, in Method of Analysis and Solutions of Crack Problems, 

G.C. Sih, ed., Noordhoff, Leyden (1973) 56-130, 


- 12 - 



5. H. Terada and M. Isida, "Analysis of Stress Intensity Factors for 
Eccentric Cracks in Plate and Strip", NAL TR-436, National Aero- 
space Lab., Tokyo, Japan (Dec. 1975). 

6. O.L. Bowie, in Methods of Analysis and Solutions of Crack Problems, 
ed. G.C. Sih, Noordhoff, Leyden (1973) 1-55. 

7. B. Gross and J.E. Srawley, "Stress Intensity Factors for Three 
Point Bend Specimens by Boundary Collocation", NASA-TN D3092 
(Dec. 1965). 

8. J.E. Srawley, International Journal of Fracture 12 (1976) RCR 475-6, 

9. J.C. Newman, Jr., "An Improved Method of Collocation for the Stress 
Analysis of Cracked Plates with Various Shaped Boundaries", 
NASA-TN-D6376 (1971). 

10. H.F. Nied and F. Erdogan, International Journal of Solids and Struc- 
tures, 15 (1979) 951-965. 

11. A.C. Kaya and F. Erdogan, International Journal of Fracture 16 
(1980) 171-190. 

12. F. Erdogan and H. Terada, "Wedge Loading of a Semi-Infinite Strip 
with an Edge Crack", Int. Journal of Fracture, 14 (1978), 399-415. 

13. M.B. Civelek and F. Erdogan, "Crack Problems for a Rectangular Plate 
and an Infinite Strip", Int. Journal of Fracture, 19 (1982) 139-159. 

14. G.D. Gupta, "The Problem of a Finite Strip Compressed Between Two 
Rough Rigid Stamps", J. of Appl . Mech., 42 (1975) 81-87. 

15. F. Erdogan, "Mixed Boundary Value Problems in Mechanics", Mechanics 
Today , Vol . 4, S. Nemat-Nasser, ed., pp. 1-86, 1978 Pergamon Press 
Inc. 

16. V.L. Hein and F. Erdogan, "Stress Singularities in a Two-Material 
Wedge", Int. J. of Fracture Mechanics, 7 (1971) 317-330. 


ACKNOWLEDGEMENT 

This work was supported by NASA-Langley under the Grant NGR 39-007-011 
and by NSF under the Grant MEA-8414477 and was completed when the second 
author was an Alexander von Humboldt Senior U.S. Scientist Awardee in 
Freiburg (i. Br.) Germany. 


- 13 - 



APPENDIX A 


The Green's Functions G. - k which appear in equations (1) and (2) 


2u (t-x)[{t-x) 2 -(y-y.) 2 J 

G xxk (x * y,y k* t) = TTra - "lTt-x)' 2 +(y-y 'I 2 ! 2 1 (k=1 ’ 2; y k =s k } 


G yyk (x,y,y k ,t) = 


G xyk (X,y, 'V t) = 


Tr(l+!C) 

L(t-x) z +(y-y k ) z J z 

2u 

(t-x)[3(y-y k ) 2 +(t-x) 2 ] 

7r( 1 +k) 

L(t-x) z +(y-y k ) z J z 5 

2u 

(y-y k )[(y-y k ) 2 -(t-x) 2 ] 

it( 1+k) 

L(t-x) z +Ly-y k ) z J z 

2u 

(y-t ) [3 ( x-x 3 ) ^*- ( t-y ) 2 J 

nN L(t-y) 2 +(x-x 3 ) 2 J 2 

2u 

(y-t)[(t-y)Mx-x 3 ) 2 ] 

t r( 1 +k) 

l_(t-y) z +(x-x 3 ) z J z 5 ( 

_ 2u 

(x-x 3 )[(x-x 3 ) 2 -(t-y) 2 ] 

L(t-y) 2 +(x-x 3 ) 2 j 2 

2 U 

— rrY— T 

(y-y j )[(y-y j ) 2 +3(t-x) 2 ] 

• r/4. ..v/w 1 .. .. TTT7 > 


(Al) 


(k=l ,2; y k =s k ) 


(A2) 


( k=l ,2; y k =s k ) 


(A3) 


(x 3 =s 3 ) , (A4) 


(A5) 


' j j-3 ; ’ 
(A7) 
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2 (y-yi)[(y-y,-) 2 -(t-x) 2 ] 

G yyj (x,y,y j’ t! = ilTbT [(t-x)-+(y-y.) 2 ] 2 ’ (J=4>5; • 

J 

(A8) 

2u (t-x)[(t-x) 2 -(y-y .) 2 ] 

GxyjU.y.y-j.t) = .rd+K) L(t - x) z+ (y-y-) z J^ ’ (j=4,5; y j =s j-3 ) * 

J 

(A9) 


2u ( x - x 6^^ x_x 6^ 2 "^ t " y ^ 2 ^ 

G xx6^ x,y,x 6» t ^ = “ ^t'O+k) "L(t-y) 2 +(x-x 6 F'F * (x 6 =s 3’ (A10) 


, v 2u ( x - x 6)^ x “ x 6^ 2+3 ^ t " y ) 2 ^ 

= tt(1+k) [(t-yj'^+Tx-Xg) 2 ? ’ (x 6 ss 3^ ’ (A11) 


2u (y-t)[(t-y) 2 -(x-x 6 ) 2 ] 

G xy6^ x ’ y ’ x 6’^ = " tt(1+k) [(t-y} 2 +(x-x 6 ) 2 J 2 » ^ x 6 =s 3^ * ^ A12 ^ 


v xk ( x .y.y k >t) - - ^n+jy ^ 1°9i:(x-t) 2 t(yy k )2] - 


(k=l,2; y k =s k ) , (A13) 


, , ^+1 y-^k (x-t)(y-y k ) 

v y k (x.y.y k .t) * ^Tfjy (~r Arctan - fr -fl^y-yj^ 


(k-1 ,2; y k -s k ) , (A14) 


(x »y »x 3 »t ) - 


1 k+1 x " x 3 (t-y)(x-x 3 ) 

tt(1+k) Arctan (t-y) z +(x-x 3 ) z} ’ 

( x 3 =s 3 ) ’ 


(AT 5) 
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V y3 (x,y,x 3 .t) - 1 ogE (t-y ) 2 +(x-x 3 ) 2 ] - u . y j^! x -pr >. 


(X3-S3) 


(AT 6) 


1 ,1+* ^3 . (x-t)(y-y 1 ) 


V x ’ y ’ y r t) = - { — Arctan ** + 


(j-4,5; yj=Sj_ 3 ), (A17) 


V yj ( x,y,y j , t) = {Sjl log[(x-t) 2 +(y- yj ) 2 ] + ( x -[)^(y ) 2> ’ 

vJ 

(j=4,5; yj=Sj_ 3 ) , (A18) 


V x6 (x,y,x 6 ,tJ ” " tt(1+ic) 1 ° 9[('t-y) 2 +(x-x 6 ) 2 ] 


+ (t-y) 2 +fx-x 6 ) 2 } * (x 6 =s 3* * 


( AT 9 ) 


i x-x fi (t-y)(x-x fi ) 

V y6 {Xsy,x 6 jt) = TTf+^T { T~ Arctan ~W + Tt^y) 2 +(x-x 6 ) 2h 

( x 6=s 3 ) . 


(A20) 
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APPENDIX B 


The kernels which appear in equation (11) 


r - u 
u xxk 7 t( 1 +k) 


^ rn xxk (x,t,c ‘^" m xxk^ H " x ’ H " t,a):]coSa ^ y " y k^ da) 

) 

(k=l ,2; y k =s k ) , 


'xx 3 tt( 1 +ic) j 


C • 
xxj 


TH+kJ I ^ m xxo^ x,t ’ a ^xxj^ H ' x ’ H " t,a ^ sina ^ y ' y j^ da * 

o 

(j-4,5; yj-Sj.j) 


'XX 6 tt( 1 +k) J 


^ m xx 6 ^ x ’ x 3 ,a ^" n, xx 6 ^ H " x ’ H “ x 3 ’ a ^ COSa ^ t " y ^ da! 

> 

(x 3 =s 3 ) , 


C yyk = JTi^T J l - m yyk (,< ’ t ’ t ‘ ) " l "yyk <H " x ’ H " t ’“ ) -' C0Se '^ y " y l< >da ’ 
0 

(k=l ,2; y k =s k ) , 


(Bl) 


^ m xx 3 ( x ’ x 3 ’ a ^xx 3 ( H ” x ’ H " x 3 ’ a ^ sinc ^ t " y ^ da ’ 

(x 3 -s 3 ) , (B2) 


(B3) 


(B4) 


(B5) 


C yy3 = [m yy3 (x,x 3 ,a)+m yy3 (H-x,H-x 3 ,a)]sina(t-y)d a , 

(x 3 =s 3 ) , (B6) 


C. • = — *■ 


J yyj tt( i+k) J 


^ m yyj^ x,t ’ a ^yyj^ H " x ’ H " t ’ ct ^ sina ^ y " y ^ da 5 



C yy6 = ^fey| ^ yy 6^’ x 3> a )- m yy6 (H ' X)H - x 3’ a):ic0Sa(t “ y)dci » 

o 

( x 3 =s 3) * (B8) 


C xyk = j Kyk (x ’ t -“ )+l "xyk (H - x ’ H - t - 0l):isl ' ncl( ^k )dc ‘ 

O 

(k=1,2; y k =s k ) , 


(B9) 


c xy3 ■ Trfcr J [i"xy3 (x ’ x 3’ !,) - m xy3 (H ‘ x ’ H - x 3-“ );lcosc ‘ (t - y)d “ ’ 

0 

(x 3 =s 3 ) , (BIO) 

CO 

c xyj Kyj (x ’ t ’“ ) - m xyj (H ' x ’ H ' t ’“ )]cos “ (y ‘ y j )d “ ’ 

O 

(j=4,5; yj=Sj_ 3 ) , (BIT) 

co 

r 

C xy6 = irXT+Ky j t m xy 6 ^ x ’ x 3’“^ 4m xy 6 < H” x ’ H ‘ x 3’“)] s i n “( t ”y) d “ • 

0 

(x 3 =s 3 ) , (B12) 


m xxk = ^ ^ (ct){[K:(e^ a ^-l ) (-l+tc 2 -2aX+2at+4a 2 tx)+4icaH(aX+at-aH)3e a ^ +x ) 

+[2ic 2 (e"^ aH -l ) (aX-at)+2aH(l-K 2 -2aH+2at+2aX-4a 2 xH+4a 2 tx)]e a ^ t ” X H, 

(k=l ,2) , (B13) 
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m xx3 = ^ ^ ( a ){E K ( e ^ a ^-l )(-l-K 2 -2aX-2aX 3 -4a 2 xx 3 )+4KciH(-l-aX+aH-ax 3 )3e a ^ x+x 3^ 

+[2tc 2 (e 2a ^-l ) 0+aX-ctX 3 )+2aH(-l-K 2 -2ctX-2aH+2aX 3 -4a 2 xH+4cx 2 xx 3 )]e a ^ x x 3^}, 

(B14) 

m .=-D ^ (a){[i<:(e 2a ^-l ) (-l-K 2 -2aX-2at-4a 2 tx)+4tcaH(-l-aX-at+aH)]e a ^ +x ^ 

+[ 2k 2 ( e 2a ^-l ) ( 1 +aX-at )+2aH(- 1 -K 2 -2aH+2at-2aX-4a 2 xH+4a 2 tx ) ]e a ^ X ^ } , 

(j=4,5) (B15) 

m xx6 = ^ ^ (ct){[<( e ^ a ^~l ) 0-K 2 +2aX-2aX 3 -4a 2 XX 3 )+4KaH(-aX-ciX 3 +aH)3e a ^ X+X 3^ 

+[2tc 2 (e 2a ^-l ) (-aX+aX 3 )+2aH(-l+ic 2 -2aX+2aH-2aX 3 +4a 2 xH-4a 2 XX 3 )3e a ^ X X 3^}> 

(B16) 

m yyk = D ‘ 1 (a){[ic( e2aH -l)(- 3 -K 2 + 6 a t +2a x - 4 a 2xt )+4 > ca H (2-a x -a t +a H )] e " a ^ t+X ^ 

+[2x 2 (e 2a ^-l ) (-aX+at+2)-2aH(-K 2 -3+2oX“6at+CaH+4a 2 Xt-4a 2 xH)]e a ^ X ^ } , 

(k=l,2) , (B17) 

niyy^D ^ (a){[i<c(e^ a ^-l )(-3+K 2 +2aX-6aX 3 +4a 2 XX 3 )+4KaH(-1+aX-aH+aX 3 )]e a ^ x+x 3^ 

+[2K 2 (e" 2aH -l)(l- a x+aX3)+2oH(-3+K 2 +2aX+6ax 3 -6aH+4a 2 Hx-4a 2 xx 3 )]e" a(x “ x 3 ) } 

(B18) 
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m -=-D ^ (a){[4e 2aH -l ) (-3+K 2 -6at+2aX+4a 2 Xt)+4KaH(-l+aX+at-aH)]e a (^ +x ) 
y y \j 

+[2k 2 ( e' 2aH -l ) ( 1 -aX+at )+2aH ( -3+< 2 +2aX+6at-6aH-4a 2 Xt+4a 2 xH ) ]e a ^ t_X ^ , 

(j=4,5) , (B19) 

m yy6 = D ^ (° t )'CC |C ( e ^ a ^“l ) (K 2 +3-2ax-6aX 3 +4a 2 xx 3 )+4KaH(-2+aX-kxx 3 -aH)]e a ^ x+x 3^ 

+[2ic 2 (e 2a ^-l ) (-2+aX-aX 3 )+2aH(-3-K 2 +2aX-6aX 3 +6aH-4a 2 xH 
+4a 2 xx 3 )]e" ct(x " x 3 ) } , (B20) 

m xyk = ^ ^ (c t ){[<(e 2a ^-l ) (T+>c 2 -2aX-2at+4a 2 tx)+4KaH(-l+aX+at-aH)3e a ^ +x ) 
+[2ic 2 (e 2a ^-l ) (-l+aX-at)+2aH(-1 -K 2 +2aH-2at+2aX-4a 2 xH+4a 2 tx)]e a ^ X ^ } , 

(k=l ,2) , (B21 ) 

^ (a){[K(e 2aH -l ) ( l-K 2 -2aX+2aX 3 -4a 2 x 3 x)+4KciH(-ax+ciH-ctX2) ]e a ( x ^ x 3) 

+[2ic 2 (e 2ct ^-l ) (ax-ax 3 ) + 2aH(l-K 2 -2ax-2ax 3 +2aH+4ct 2 x 3 x-4a 2 Hx)3e a ^ x x 3^}> 

(B22) 

m .=D ^ (a){[ic(e 2a ^-l ) (-k 2 +1 -2aX+2at-4a 2 tx)+4icciH(-aX-at+aH)]e a ^ +x ^ 
xyj 

+[2K 2 (e" 2 “ H -1)(ax- a t)+2aH(l-.c 2 +2aH-2c l t-2ax-4a 2 xH+4a 2 tx)]e a(t ' x) }, 

(j=l,2) , (B23) 
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m x yg=-D ^ (a){[K(e 2 ^-l ) (-l-K 2 +2aX+2aX2-4a 2 XX2)+4icaH(l-aX-aX2+aH)Je a ^ x+x 3^ 
+[2K 2 (e" 2aX -l ) (l-aX+c 1 X3)+2c i H(l+ K 2 -2aH+2aX3-2aX+4a 2 xH-4c t 2 XX3)]e" a(x ‘ X3j } , 

(B24) 

D(a) = ic 2 [e 2aH +e"*^ aH -2]-4a 2 H 2 . (B25) 
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Fig. 1 The basic crack geometry considered. 




Fig. 2 Normal stress along the clamped boundary in a plate without 
a crack (k=1 .8) . 
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Fig. 4 Normalized Mode II stress intensity factor at the comers of 
a clamped rectangular plate without a crack (k=1.8). 
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Fig. 5 Normal stress along the clamped boundary in a plate with a 
central crack (H/2s-j=l), <=2.2. 
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